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We derive a phase diagram for the pseudogap onset temperature T* (associated with the breakdown 
of the Fermi liquid state, due to strong pairing correlations) and the superconducting instability, 
T c , as a function of variable pairing strength. Our diagrammatic approach to the BCS - Bose- 
Einstein cross-over problem self consistently treats the coupling between the single particle and pair 
propagators, and leads to a continuous evolution of these propagators into the standard T < T c 
counterparts. A rich structure is found in T c which reflects the way in which the superconducting 
instability at T c is affected by the pseudogap A pg . An important consequence of Cooper-pair- 
induced pseudogaps is that the magnitude of T c is sustained, even when A pg > T c . 
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It is generally agreed that the pseudogap state of 
the underdoped cuprates represents some type of pair- 
ing above T c which is postulated to derive from neutral 
spinon pairs (lj, spin or charge jl| density wave states, 
or from some form of (2e) Cooper pairing which foreshad- 
ows the ultimate superconducting state. Each of these 
scenarios must not only address the nature of the exotic 
(ie., pseudogapped) normal state but also the transition 
from this state to its associated superconducting instabil- 
ity. Indeed, there is an extensive literature on the prob- 
lem of establishing superconductivity in the presence of 
gaps in the electronic spectral function (J]. In this pa- 
per we investigate this issue under the presumption that 
the normal state spectral function gap is associated with 
Cooper pairing above T c . Moreover, by self-consistently 
coupling the single particle and pair properties we es- 
tablish the nature and physics of the pseudogap state. 
In this coupled scheme, we demonstrate how gapping in 
the electronic spectrum affects T c and argue that this 2e 
pairing is not particularly deleterious to the superconduc- 
tivity, despite the relatively large size of the pseudogap 
(Apg > T c ). Our results are consolidated into a phase 
diagram in which the pseudogap phase occupies a large 
temperature range in the limit of moderately strong su- 
perconducting coupling. 

An important contribution of this paper is to revisit 
the BCS - Bose Einstein cross-over problem || using a 
conserving diagrammatic "mode coupling" formulation. 
Above T c this approach leads to one and two electron 
Green's functions which are formally continuous with 
their below T c counterparts obtained in the standard the- 
ory of superconductivity. This approach also reproduces 
the expected limits associated with the extreme weak and 
strong coupling regimes of previous saddle point schemes 
|^| . The intermediate coupling regime is of particular 
importance to the cuprates and Monte Carlo simulations 
suggest that deviations from Fermi liquid behavior 
are present there. Our previous diagrammatic calcula- 
tions B indicate that long lived pairs introduce a gap 



in the electronic spectral function by blocking available 
single particle states around the Fermi surface. This 
gap correlates with a regime of resonant pair scattering 
which forms a natural intermediate state between the free 
fcrmions of the weak coupling, and the composite bosons 
of the strong coupling limits. In the present fully self 
consistent theory, this resonant scattering will, itself, be 
modified by pseudogap effects. 

We consider a generic Hamiltonian consisting of 
fermions in the presence of an attractive interaction, 
which may be modeled as Vk,k' = g^k^k 1 , where tp^ = 
(1 + k 2 /fco)~ 1//2 an d 9 < is the coupling strength of- 
ten expressed in units of g c = —Air/mko |9J. The BCS - 
Bose Einstein cross-over problem is characterized by cou- 
pled equations for T c and the chemical potential, \x. The 
former derives from the Thouless criterion for the pair 
propagator or T-matrix, and the latter from the parti- 
cle number equation, which involves the single particle 
propagator. In general, the full treatment of the cou- 
pling between the pair and single particle propagators is 
referred to as "mode coupling" . In the present formula- 
tion the coupled equations are obtained following earlier 
literature on superconducting fluctuation effects [Efj 
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where the Green's function is given by G, ■ - w, 

S k,^,' G k^ = iu >l ~ £ k; are the even/odd 

Matsubara frequencies, and the electronic dispersion is 
£k = k 2 /2m — fi. 

At high temperatures where mode coupling or "feed- 
back effects" may be ignored || , all Green's functions in 
Eqs. [I] and ^| may be replaced by bare propagators. If we 
restrict our attention to moderately strong coupling con- 
stants, g/g c ~ 1 (associated with fcp^/27r ~ 1 — 10), we 
find a Fermi liquid phase above a crossover temperature 
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T* . The break-down of the Fermi liquid at T Res ~ T* co- 
incides with the onset of resonant pair scattering. As the 
temperature is lowered below T* the gap in the spectral 
function grows as the resonant states become more long 
lived. Finally, at T c , coherent pairs condense into a su- 
perconducting ground state. The value of T c is obtained 
by solving Eqs. [I] and |2|, which, at T c , become 
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where £k — 



The pseudogap parameter 



is, at arbitrary temperatures, defined to be 



A 2 

pg 
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(5) 



where /(w), b(uS) = (exp(w/T) ± 1) 1 . It is evident that 
A 2 g in Eq. || coincides with the square amplitude of pair- 
ing fluctuations, g 2 < Sacc >. 

It should be noted that Eqs. || and [| provide an im- 
portant validation for the diagrammatic scheme of Eqs. |l| 
and |^. These two equations show that a pseudogap in 
the normal state spectrum preserves the overall "BCS- 
like" structure of the T c equation as well as the diagonal 
component of the single particle Green's functions. This 
continuity implies that, even though these equations are 
derived above T c , they are formally identical to their one 
and two particle counterparts of the standard supercon- 
ducting state jTTJ . This is in contrast to the use of fully 
rcnormalized Green's functions everywhere. 

In order to evaluate T c it is necessary to further sim- 
plify Eq. ||. As a first step, we investigate the T-matrix 
in the presence of a pseudogap. To capture the essen- 
tial physics, we use a time dependent Ginzburg-Landau 
(TDGL) formulation 
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where the TDGL parameters can be calculated by a for- 
mal expansion of the full T-matrix at low frequencies and 
momenta. For intermediate coupling and at sufficiently 
low n and q, the T-matrix assumes the resonant form 



f2 — fL 



(7) 



where ao is predominantly real. It follows from Eq. |^, 
that feedback effects, associated with a gap in the single 
particle spectrum, lead to a gap in the imaginary part of 
the inverse T-matrix at T c . This latter gap affects T q in 



an important way: when il q < 



A r 



0. This gap 



is slightly smeared out for temperatures between T* and 
T c , but, nevertheless T q remains small. Thus mode cou- 
pling leads to a stabilization of resonance effects, and, 



thereby an amplification of pseudogap behavior. This 
TDGL formulation is in contrast to results obtained from 
previous cross-over TDGL theories ]l2"| which omit mode 
coupling, as well as from applications to dirty supercon- 
ductors JTJJ], where mode coupling is included, but gap 
effects arc absent. 

The above results can now be applied to rewrite Eq. |5| 
and, thereby, compute the superconducting instability 
temperature, T c . The value of T c is obtained from a so- 
lution of the Thouless condition (Eq. and the number 
equation in the presence of the self-consistently deter- 
mined pseudogap amplitude A pg . The coupled equations 
for this transition temperature follow from Eqs. |^-|| and 
are given by 
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where v 2 . = (1 



e k /£k)/2. 

We solved Eqs. [s||lO| numerically to obtain T c as a func- 
tion of g I g c , along with the resonance onset temperature 
TRes, which we associate with T* . Our results are plotted 
in Fig. |l|. Also shown (dashed line) is the superconduct- 
ing instability temperature T computed in the absence of 
mode coupling. When there is an appreciable pseudogap, 
T c and T* vary in an inverse fashion with g/g c - This is a 
consequence of the pseudogap which suppresses T c , but 
is not present at T*. In this way the pseudogap regime 
is greatly enhanced by mode coupling contributions. 

The overall behavior of T c is compared with that ob- 
tained from the saddle point approximation, as well as 
that of strict BCS theory in Fig. ||(a). Here the inset 
plots the chemical potential and pseudogap function (at 
T c ), along with T c . It can be seen from the inset, that 
the maximum in T c is associated with /i ~ A pg and the 
minimum with ji = 0. Indeed, the complex behavior of 
T c shown in Fig. || can be understood on general physi- 
cal grounds. A local maximum appears in the T c curve 
as a consequence of a growing (with increased coupling) 
pseudogap A pg in the fermionic spectrum which weak- 
ens the superconductivity. However, even as A pg grows, 
superconductivity is sustained. In the present scenario, 
superconductivity is preserved by the conversion of an 
increasing fraction of fermions to bosonic states, which 
can then Bose condense. Once the fermionic conversion 
is complete (/x = 0), T c begins to increase again with cou- 
pling. This is a consequence of the decreasing pair size 
and concomitant reduction of the Pauli principle repul- 
sion. 

The behavior of T c on an expanded coupling constant 
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scale, for different ranges of the interaction (parameter- 
ized by ko/k-p) is shown in Fig. |^(b). The limiting value 
of T c for large values of g/g c approaches the ideal Bose- 
Einstein condensation temperature Tbe = 0.218-Ep as 
ko — ► oo Jl3|. The qualitative shape of the T c curve, how- 
ever, is retained as long as k /kp is greater than about 
0.5. For larger range interactions, the solution disap- 
pears for some g/g c ] then, when jj, is sufficiently nega- 
tive, the transition reappears approaching a continuously 
(with fco/fcp — » 0) decreasing asymptote. 

In order to implement the mode coupling scheme away 
from T c , we introduce a two parameter fit to the full 
self-energy which is based on, and evolves smoothly into 
Eq. I at T c , 



rime 



A p>£ 

+ £k + H ' 



(11) 



This model makes it possible to do systematic numerics 
in the well established pseudogap phase, but it is not 
appropriate in a regime in which the pseudogap breaks 
down as it evolves towards the Fermi liquid state. The 
instability of the Fermi liquid, or pseudogap onset T* , 
is, however, addressed using the lowest order expansion 
of Eqs. lland ^, as discussed above. The parameter A pg 
in Eq. IT] is defined from Eq. |j| which is numerically 
implemented via 



— Im S L 



(12) 



Similarly, the parameter 7 is obtained from the height of 
the peak in —Im Ek, w which is defined to be A^ g ip^/j. It 
should be stressed that at T c the parameter 7 in Eq. |ll| 
becomes zero, as is implied by Eq. |3|. This results from 
the continuity of the pseudogap and superconducting 
phases, and can be seen explicitly by noting that T c oc- 
curs when the T-matrix becomes singular at 57 = and 
q = (the Thouless criterion). As a result, a divergent 
peak in Im £k,u occurs at u> = — e^. 

The coupled equations for Ek,u an d ^ q q? as a func- 
tion of temperature, were solved numerically using the 
model self-energy of Eq. [ll] and the TDGL approxima- 
tion of Eq. ||. In a given iteration A pg and 7 were used to 
compute tq t Q via Eq. ^ which was then approximated by 
Eq. H for use in Eq. |]. The output of the latter was then 
used to extract A pg and 7 via Eq. ^ and the procedure 
repeated until convergence. In this way A pg , 7, f2 q , and 
T q were all determined self-consistently as a function of 
temperature and coupling constant. 

To emphasize the important coupling between these 
electronic and pair properties, we compare the evolution 
of the single particle parameters A pg and 7 with the pair 
parameters Sl q= o and r q=0 m Fig- R f° r 9 1 9c — 1-2. The 
associated spectral functions over this range of tempera- 
tures are shown by the inset in Fig. [|(a). The shaded re- 
gions indicate where the two parameter model self-energy 



of Eq. [Tl] breaks down, as does the lowest order theory; 
even though the pseudogap persists, the system begins 
to cross-over towards the Fermi liquid. Thus the param- 
eters A pg and 7 are indicated only for a limited range of 
temperatures. The temperature for Fermi liquid onset, 
Tros ~ T* , is just beyond the shaded region corre- 
sponding to the regime of validity of the lowest order 
theory. This temperature can be read off from the end- 
point of the two curves in Fig. |](b). Moreover, the ratio 
Ap g /7 which represents the peak height of /mS^u, ap- 
pears to extrapolate smoothly from the mode coupling 
results to zero at T^ es , at which point the Fermi liquid 
state sets in and Im £k,w changes character. As indi- 
cated in Fig. ||(a), the disappearance of the pseudogap 
arises via a reduction in A pg , while the peak broadening 
remains relatively small. This latter is a consequence of 
the extended lifetime of the pairs due to a full (s-wave) 
gap in the T-matrix. This gap is also responsible for 
the relatively small values (when compared to conven- 
tional TDGL theory jio]]) of r q= o, which is associated 
with the inverse pair lifetime, shown in Fig. ||(b). At 
higher temperatures r q= o and f2 q= o are plotted to join, 
after extrapolation, onto the results of the lowest order 
theory. 

In summary, within a BCS - Bose Einstein cross-over 
picture, we have presented a quantitative phase diagram 
which compares the temperature onset of the pseudogap 
with the onset of a novel superconducting state, associ- 
ated with pseudogapped fermions. Mode coupling effects, 
which were important for this analysis, considerably en- 
hance the pseudogap regime. We have demonstrated (for 
the s-wave case) that this pseudogap disappears with 
temperature, as it evolves into a Fermi liquid state, prin- 
cipally by a reduction in the gap size A pg . While we have 
not established detailed connections to the cuprates, the 
inverse of g/g c in Fig |l| can be loosely associated with 
hole concentration jwj. Of importance for the cuprates 
is the prediction that the effective inverse lifetime in the 
electronic self-energy 7, which can be deduced experi- 
mentally jig] , varies continuously to zero at T c . Finally 
changes in our results associated with the d-wave, lay- 
ered structure of the high T c systems can be anticipated: 
the nodal structure of d-wave pairs will weaken the gap 
in the T-matrix which played an important role in de- 
termining the detailed temperature evolution of A pg (T) 
and 7(T) above T c . Moreover, quasi-two dimensional- 
ity will considerably lower the energy scales and enhance 
the pseudogap regime, particularly as the insulator is ap- 
proached. Despite these omissions, our physical picture 
of the interplay of the pseudogap and superconducting 
instability is expected to be qualitatively general, within 
a precursor superconductivity scenario. 
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FIG. 1. Pseudogap phase diagram, which indicates the 
Fermi liquid breakdown co-rresponding to pair resonance on- 
set, at Tr,os, the fully self-consistent T c and the transition 
temperature T in the absence of mode coupling. 
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FIG. 2. (a) Comparison of the T c curves in BCS theory, 
the Gaussian approximation of Nozieres and Schmitt-Rink 
(NSR), and our mode coupling theory (MC). The inset plots 
the corresponding values of fi, A pg , at T c . (b) Variation of T c 
for three different pairing interaction ranges. 
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FIG. 3. Evolution of the parameters which characterize the 
electronic self-energy of Eq. llL (a) , and the pairing resonance 
(b). Shaded region represents the breakdown of Eq. [ll] as well 
as lowest order theory. The latter is valid to right of shaded 
region. The values of Ap g /7 were divided by a factor of 10 
and g/g c = 1-2. Spectral functions are plotted in the inset, 
for the three temperature regimes. 
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